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Abstract 

Electronic transport in a triple quantum dot shuttle device in the presence of an ac field is 
analyzed within a fully quantum mechanical framework. A generalized density matrix formalism 
is used to describe the time evolution for electronic state occupations in a dissipative phonon bath. 
In the presence of an ac gate voltage, the electronic states are dressed by photons and the interplay 
between photon and vibrational sidebands produces current characteristics that obey selection 
rules. Varying the ac parameters allows to tune the tunneling current features. In particular, we 
show that coherent destruction of tunneling can be achieved in our device. 

PACS numbers: 
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The study of nanoelectromechanical systems (NEMS) has drawn a great deal of interest, 
both at an applied and at a fundamental level, due to novel transport regimes arising from the 
strong interplay between electrical and mechanical degrees of freedom 1 . This is of particular 
interest in the field of molecular electronics where, for example, the vibrational modes of 
a Cqo molecule are coupled to its charge state giving rise to new features in the current 
through the molecule 2 . 

Shuttling is an example of nonlinear transport regime in NEMS, where an oscillating 
nano-grain acts as a shuttle carrying one electron from the left to the right contact. A 
paradigmatic proposal of such a device was made in the seminal work by Gorelik et al. 3 , 
and, since then, there have been different models that describe the mechanical degree of 
freedom using semiclassical 4 or quantum mechanical 5,6 approaches. In the presence of an 
ac electric field, it is known that the interplay between the ac and the natural frequency of 
the device gives rise to rich dynamics with strong nonlinear characteristics. These can be 
modified by tuning the ratio of the two frequencies or the intensity of the applied field. Up 
to now, this problem has been addressed for the case of an oscillating quantum dot (QD) 
using a semiclassical approach to describe the oscillations 4 . 

In this letter we introduce a full quantum mechanical description to study the dynamics of 
a triple dot quantum shuttle (TDQS) with an ac gate voltage applied between the different 
QDs and including quantum optical damping as a dissipation bath 5 . Using the master 
equation approach for the reduced density matrix (RDM), we show how tunneling current 
resonances occur through different combinations of photonic and phononic sidebands and 
derive sum rules to explain these resonances. These features strongly depend on the ratio 
between the frequency and intensity of the ac field. Thus, by tuning the ac frequency we 
can obtain information on the vibrational frequency of the resonator. Finally, we find that 
the ac parameters can be manipulated so that coherent destruction of tunneling 7 ' 8 (CDT) 
can be observed. 

We consider a TDQS, which consists of an array of three QDs where the central movable 
dot is flanked by two static dots at fixed positions ±x . The oscillation of the central dot 
will affect the tunneling rates, which are now position dependent. The charge transport is 
within the Coulomb blockade regime, so the charging energy is assumed to be sufficiently 
large that only one transport electron can occupy the chain of three dots at any given time. 
Our model introduces a time- dependent ac potential V(t) = V ac cos(u ac t), applied between 
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the left and right QDs, where V ac and uj ac are the amplitude and frequency of the applied field, 
respectively (see Fig. 1). The system is modeled by means of a Hubbard like Hamiltonian, 
H = H TQS + H leads, where the unperturbed hamiltonian is H TQS = H + H osc + H tun + H ac , 
where H = XL £ s c I c s represents the Hamiltonian of the left-hand (/), central (c), and 
right-hand (r) dots, for spinless electrons with respective localized states e s (s = l,c,r). 
The energy of the central QD is e c (x) = [si — (eb/2xo)(x + xq)], where it is assumed that 
it undergoes a Stark shift proportional to its position due to the voltage bias across the 
device e b = (e t — e r ); H osc = hwSd is the energy associated to the oscillation of the central 
QD, with frequency uj. H tun = ^2 s '=i r ^s'(^) \ c \ c s' + h.c], represents the tunneling between 
all three dots, with position dependent tunneling rates T] = —V exp~ a ^ Xo+x \ and T r = 
—V exp~ a ( x °~ x \ where V is the tunneling amplitude and a is the inverse of the tunneling 
length. The position operator x measures the displacement of the vibrational mode and is 
given by x = Ax zp + d), where Ax zp = (h/2mu) 1 ^ 2 is the zero-point uncertainty position 
of the oscillator. H ac = (V ac /2) cos(uj ac t) (c|q — cjc r ), provides the effect of the ac potential, 
which is introduced as an oscillation of opposite phase on the energy levels of the right 
and left dots flanking the central dot and Hi ea d s = J2 a ^2k a a- e k a c k a a c k a v The equation 
of motion of the system using the standard RDM approach 9 , generalized to include the 
environment of the oscillator, is given by: p^ s = (p 1 J s )rdm + {p l j s )diss-> where p^ s represents 
the matrix element with the dot states q and s (q,s — 0, 1, c, r), with the oscillator states i, 
and j = 0, 1,2, ...). The term (p^ s )rdm incorporates transitions between the leads and 
the outer dots: 

(p%)rdm = -jjfi, p}% 

+ j Ed^ q ( T d^ q Pdd - r ?-d p&);q = s 
\~i (E^ 9 r 9 _ d + J2 d ^s T s^ d ) p%; q^s. 

Here we have also considered the transitions r r _^ = r ^ = T. The term [p^ s ]diss accounts 
for the dissipative effects of the oscillator's environment. We assume a damping model 
composed of a bath of oscillators at a fixed temperature T, to which the vibrational mode 
is weakly coupled 10 . In the limit of small values of T, the matrix elements are given by 
[p%}diss = + j)/2)]p% + 7[(i + l)(j + l)] 1/2 p$ 1J+1 , where 7 is the classical damping 

rate of the oscillator. Finally, the matrix elements of the tunneling amplitudes T\ 3 r in the 
oscillator basis can be computed. 
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In what follows, we shall explore the main features of the current I ac for a TDQS driven 
by an ac gate voltage. We have considered weak coupling between the three QDs, a regime 
characterized by large values of a. The calculations are performed by numerical integra- 
tion of the RDM and averaging over the resulting electronic current measured at the right 
QD, i.e., I ac = eT J2 i \p™ r ] av . We begin our analysis by exploring I ac as a function of the 
bias Eb for different values of the ac parameters. All of the calculations are worked out in 
units where h = 2m = uo = e = 1 and N = 6 oscillator states were included. For the 
particular values of the voltage bias across the device £& considered in our calculations, we 
have numerically verified that including more oscillator states does not significantly change 
the current characteristics for typical experimental parameters. In Fig. 2(a) we plot the 
time-averaged current I ac (solid line) for a particular value of the ratio £ = (V ac /2uj ac ) = 3.0, 
where ui = pu) ac (p integer). We can appreciate a much richer structure for I ac than the 
undriven case (dashed line). We find a series of peaks or maxima appearing at specific values 
of Eb, which are integer multiples of u ac , and obeying the following sum rules: 

Eb + vhw ac = nftw, (2) 

(£&/2) + v?Vjj ac = nhw. (3) 

Here, n corresponds to the vibrational level at which the transmitted electron ends up in 
the right dot, and v to the number of emitted or absorbed photons needed to reach that 
vibrational level. The peaks that satisfy the sum rule of Eq. (2) are positioned at a value 
of Eb equal to an odd multiple of u ac , and are identified by [v,n]. These correspond to an 
alignment between photon-assisted tunneling (PAT) sidebands of the vibrational resonator 
states of the right and left dots. As in the undriven case, these peaks are suppressed in the 
very weak coupling regime. The peaks satisfying Eq. (3), positioned at an even multiple of 
u ac , are identified with (v, n), and correspond to alignment of sidebands of all three dots. In 
order to further corroborate the label assignments, we calculate the occupation probabilities 
for the right QD p™ r , which determine I ac . In Fig. 2(b), we present the results corresponding 
to the peak at E b ^ 0.4 [with labels (—1, 0), (+4, 1)], where we can clearly see that the most 
important contribution arises from Similarly, for all the resonances we have verified that 
the most important contribution comes from the satellite with the least number of photons, 
which is consistent with a low field intensity. 
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In what follows we shall introduce a useful transformation in order to demon- 
strate that the different peaks associated to the PAT current are in fact associated 
with the emission or absorption of v photons. Performing a unitary transformation 11 



U = exp 



-i(V ac /2u ac ) sin(u; a( i)(cjq - c\c r ) 



on the Hamiltonian H, leads to a trans- 



formed tunneling Hamiltonian H' tun with time dependent tunneling amplitudes T} = 
Ti E^(-1)"^(0^, and f r = T r EZ-Z (^Y MO ■ Thus, the tunnel- 

ing amplitudes become renormalized by a sum of Bessel functions of order u, where v is the 
number of absorbed or emitted photons. By computing the current using this approach we 
obtain an excellent agreement with our previous calculation. The peaks in I ac can be ex- 
plained by an interplay of the vibrational and photoassisted sidebands as follows: the TDQS 
can be seen as a system of three fixed QDs, whose energy levels are separated by (e:&/2), 
and the effect of the oscillator is to split these energies into vibrational satellite sub-bands 
separated by the vibrational energy hw. When the ac voltage is applied across the device, 
an additional level splitting of their energies occurs by an amount of friv ac , giving rise to new 
conduction channels (resonances) which are accurately described by the sum rules. 

We can now exploit the analytical properties of the Bessel functions in order to control 
the current characteristics by an adequate manipulation of the ac parameters. In Fig. 3 we 
use £ = 2.405, for oo ac = 0.2 and a lower V ac = 0.96192, corresponding to the first zero of 
Jo(0- I 11 this case, direct tunneling processes involving v = at e& ~ 0.0, 2.0, are suppressed 
due to CDT 7 ' 8 . Here, the energy levels of the vibrational subbands of the three dots are 
aligned, but 7} ks and T r ps 0. Note that these values of Sb also correspond to resonances 
in the undriven case. Interestingly, the current peaks at Sb ^ 0.6, 1.0, and 1.4 [with labels 
[v, n]], observed in Fig. 2, are also significantly diminished. This decrease in the current is 
due to two factors: the small overlap of the wave functions in both left and right QDs and 
the decrease in the intensity of the ac field, which yield a small probability for processes 
involving several photons. We see that the main peaks are governed by the emission and 
absorption of several photons [y = ±4, ±1) and correspond to alignment of the levels of 
all three dots. 

In conclusion, we have analyzed the characteristics of the electronic current in an ac 
driven TDQS device using a full quantum mechanical approach. Our results show that 
by tuning the ac field parameters we can control the current contribution from differently 
vibrational modes. This constitutes a new mechanism to obtain information on vibrational 
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modes through current measurements as for example in the context of molecular shuttle 
devices. 
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FIG. 1: Driven TDQS where oscillating tunneling barriers between the central quantum dot and 
its neighbors are considered. The system is driven by an ac potential V(t) = V ac cos(uj ac t) and the 
left and right tunneling barriers are rigidly coupled to the leads. 

FIG. 2: (Color online) (a) I ac (in units of V) vs e (blue solid line) for the ac parameters u) ac = 0.2, 
and V ac = 1-2, and typical values of the parameters {V, 7, xq, T, a}={0.5, 0.01, 5.0, 0.05, 0.4}. 
The red dashed line corresponds to I ac in the undriven case. The ac current peaks are labeled with 
[u,n], and (i/,n), associated to the sum rules given by Eqs. (2), and (3), respectively. In (b) we 
show the occupation probabilities of the right dot p** r , for i = (green line), i = 1 (red dashed 
line), and i = 2 (blue dotted line), as a function of time, corresponding to the current peak at 
e b ~ 0.4. 

FIG. 3: (Color online) I ac (in units of Y) as a function of (blue solid line), for uj ac = 0.2, 
V ac = 0.96192 (£ = 2.405), corresponding to the first zero of Jo(£)j showing CDT. Parameters for 
the undriven case (dotted line) are the same as in the previous figure. 
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